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Abstract. This paper presents results from a flow visualization study of the 
wake structures behind solid spheres rising or falling freely in liquids under the 
action of gravity. These show remarkable differences to the wake structures 
observed behind spheres held fixed. The two parameters controlling the rise or 
fall velocity (i.e., the Reynolds number) are the density ratio between sphere and 
liquid and the Galileo number. 



1. Introduction 

This year's cover illustration shows the boundary layer separation from a light sphere 
rising in quiescent water at high Reynolds number. The wake structure behind the 
sphere is visualized with the Schlieren technique. The experiment is part of a project 
to study wake structures behind falling and rising spherical particles in a large range 
of Reynolds numbers. 

In the past years extensive numerical investigations [1-5] have established how the 
wake of a sphere held fixed in a uniform flow undergoes a series of transitions as the 
Reynolds number Re — U d/v is increased. Here U is the free stream velocity, d the di- 
ameter of the sphere, and v the kinematic viscosity of the water. It was found that the 
wake is axially symmetric up to Re = 212. Above this value a planar-symmetric wake 
is found that consists of two steady counter-rotating threads. At Re ~ 270 there is a 
further transition and the planar-symmetric flow becomes time-dependent: opposite- 
signed streamwise vortices then form a series of loops that resemble hairpin vortices. 
As the Reynolds number is further increased, the flow gradually becomes more irreg- 
ular and finally turbulent. The Digital Particle Image Velocimetry (DPIV) measure- 
ments by Briicker [6] and the flow visualization studies by Schouveiler & Provansal 
[7] have confirmed most of these numerical results and have further elucidated the 
sequence of transitions. 

For freely moving spheres the Reynolds number is defined by the measured mean 
velocity of rise or fall of the sphere Ut and the corresponding 'Reynolds number' 
becomes Rer = {UT)d/v. The mean velocity is the time averaged velocity of the 
sphere, not including the acceleration of the sphere from rest. The flow looses its 
axial symmetry at a critical Reynolds number which is not significantly affected by the 
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Figure 1. Top view of tlie Schlieren set-up used to visualize the sphere wakes. 



density ratio ps/p [8]: i?ecr=211.9 for ps/p oo (i.e. sphere held fixed), Recr = 206.3 
for ps/p = 0.5 and Rccr = 205.8 for ps/p = 0.0 . This is in good agreement with, for 
example, the experimental results on solid spheres [9], on surface-contaminated gas 
bubbles [10], and on wake visualizations in experiments with drops of tetrachloride 
and chlorobenzene falling in water [11 13]. As pointed out by Natarajan & Acrivos 
[14], these drops must have behaved effectively as solid spheres due to presence of 
surface-active impurities, and these visualizations have therefore often served as a 
basis of comparison with numerical studies on fixed spheres. 

What happens for freely falling or rising spheres at higher Reynolds number, 
which is more common in multiphase flow applications? How are the wake structures 
and transitions observed for the fixed sphere case modified? Is there a (clear) difference 
in wake structure between rising and falling spheres? In this paper we present flow 
visualizations of the wakes behind freely moving solid spheres at large Reynolds 
number {Re = 450 — 4623) for which the density ratio ps/p is in the range 0.50 
to 2.63. 

2. Experimental details 

The flow visualizations of the sphere wakes were carried out in a transparent tank 
(0.15 x 0.15 X 0.5 m^) filled with decalcified water. Smooth plastic spheres with 
diameters between 1.5 mm and 10 mm and densities between 500 kg/m'^ and 2781 
kg/m^ were released from rest. By means of an optical system consisting of two 
LED-lights, pinholes, lenses and mirrors, two perpendicular images of the particle and 
its wake were created and recorded at 500 frames/s with a CCD-camera (figure 1). 
Hence, each image consists of two perpendicular views of the same sphere. The images 
are taken at a position in the transparent tank where the spheres do not accelerate 
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Table 1. Parameter values in our visualizations; d in mm, ps in kg/m^. G 
is defined by equation (1) and Rbt is the mean Reynolds number {Rbt = 

(UT)d/(u)). 
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anymore. The wake was visualized using the Schheren technique. To this end a 
small vertical temperature gradient in the water was maintained (1 K/cm.). The 
mean water temperature at the measurement section was 302 K, with corresponding 
values of the density ((p)) and viscosity ((;/)) of 996 kg/m'^ and 0.802-10"^ m^/s, 
respectively. Hence, the Reynolds number is based on the mean viscosity and is defined 
as Rct = {UT)d/ {v) . It turned out to be difficult to keep a constant temperature 
gradient. Therefore the error in the mean water temperature at the measurement 
section is about 3 K, leading to a relative error in the viscosity of 10 %. As opposed 
to the fixcd-sphcrc problem, the Reynolds number for freely moving spheres is not 
an independent parameter. Following Jenny et al. [8] we choose as independent 
dimensionless variables the ratio Ps/{p) of the densities and the Galileo number 

{\{pj{p)-l)\gf'dy^ 

1/2 

Since {\{ps/{p) — can be considered as a velocity scale, G plays a similar 

dynamical role as the free-stream Reynolds number in the case of fixed sphere. The 
parameter values for which we made the flow visualizations are summarized in Table 
1. 
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Figure 2. Stereoscopic views of falling and rising spheres and their wakes. The 
left part of each frame shows the xz-plaiie and the right part the yz-plane. In 
each case the sphere diameter is 6 mm. The values of the parameters Ps/{p), G, 
and RcT are, respectively: (a) 1.04, 359, 546; (6) 1.05, 394, 608; (c) 0.96, 355, 
647; (d) 0.95, 390, 656. 



3. Observations 

Figure 2 shows stereoscopic images of the wake structure behind falhng spheres 
with densities approximately 4% (fig. 2a) and 5% (fig. 2h) higher than that of the 
surrounding liquid, and, for comparison, that behind rising spheres with densities 
that are approximately 4% (fig. 2c) and 5% (fig. 2d) lower. In all these examples the 
sphere diameter is 6 mm, so that the parameter G is roughly identical in cases (a) and 
(c), and in cases (6) and [d). The lighter spheres have a slightly higher vertical velocity 
than the heavier spheres, as indicated by Rer in table 1. The wakes of the falling 
spheres appear to have a more 'irregular' structure, and the path followed by these 
spheres shows much larger deviations from a straight vertical line. These pictures 
illustrate that the density ratio Ps/{p) matters, even at values close to 1. 

Figures 3 a and 3& give examples of a phenomenon that we believe to be 
characteristic for spheres following a zigzag path, namely that the two counter-rotating 
threads in the wake cross at the centerline of the zigzag (indicated by '1' in fig. 3 a). 
The presence of these threads of opposite-signed streamwise vorticity implies that the 
sphere experiences a lift force. As a consequence of the periodic crossing of the threads 
this force is always directed towards the zigzag center- line (see the sketch in fig.3c). 
A similar observation was made by De Vries et al. [15] on the wake behind zigzagging 
gas bubbles. 

Schouveiler & Provansal [7] remark that for a fixed sphere "the dynamics of 
the two opposite-sign streamwise vortices ... presents a striking similarity with the 
long-wavelength (or Crow) instability of a pair of counter-rotating parallel vortices" 
and further "such a vortex pair instability could be responsible of the appearance 
of unsteadiness in the sphere wake". Figure 3 suggests that the situation is slightly 
different for freely moving spheres. Here it appears that close to the sphere each of 
the vortices first develop a 'kink' (indicated by '2' in fig. 3a), a process in which the 
curvature of the vortices presumably plays an important role [16] . As the kinks develop 
further downstream of the sphere they come near each other and finally combine into 
what resembles a hairpin vortex (indicated by '3'). This sequence of events can be 
seen in the flow visualizations presented in figure 4, see also figure 6 of ref. [11]. 
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As the kinks develop and hairpin-like vortices are formed further downstream, a 
pattern results. Lunde & Perkins [17] interpreted this pattern as a scries of hairpin 
vortices of alternating sign, shed periodically by the spheres at the extremes of the 
zigzag path. Our visualizations suggest instead that the streamwise vorticity produced 
at the surface of the sphere docs not change sign; the legs of the like-signed hairpin 
vortices cross at the centerline of the zigzag. 

Figure 26 is an example in which more than one kink develops in a half-period 
of the zigzag. We have not yet been able to determine the conditions (in terms of the 
parameters Ps/{p), G, or Rer) that select the number of kinks that are formed. What 
is remarkable is that the development of the kinks and the subsequent formation of the 
hairpin vortices do not seem to affect the trajectory of the sphere. This corroborates 
the opinion that at high Reynolds numbers the details of the vorticity distribution 
very close to a body basically determine the forces that it experiences. 

We will now turn to experiments with density ratios more different from one. Recently 
Jenny et al. [18] reported on their numerical work on freely moving spheres in a 
Newtonian fluid. They focussed on the frequencies in the wake and the path of 
the sphere in the parameter space spanned up by the density ratio and the Galileo 
number. Figure 5 reproduces their phase diagram. The numbers in the diagram 
refer to the numbers of the experiments given in table 1 . A lot of our experiments are 
outside their investigated region and new experiments should be done in the interesting 
regions around a Galileo number of 200. Further experiments should focus more on 
the frequencies in the wake of the sphere and compare this to the frequencies given 
by Jenny et al. (see caption of figure 5). Furthermore, we must stress that wake 
visualizations with the Schlieren method demand a temperature gradient in the water. 
Hence the density and viscosity of the water are not constant through the entire flow 
field and the local Galileo number will not be constant. The differences between the 
mean Galileo number and the local Galileo number can reach 10 % and must be taken 
into account when analyzing figure 5. 

A striking difference between our experimental data and the numerical data of 
Jenny et al. is the behavior of falling spheres with a density ratio close to one. From 
figures 2 and 8 it can be seen that these falling spheres can also fall in a non- vertical 
path. This contradicts Jenny et al. who claim that only rising spheres can go in a 
zigzagging motion (the circles in the phase diagram figure 5). 

From our experiments one concludes that for increasing Reynolds number the 
wake becomes more irregular (figures 6 to 8). The two-threaded wake structure is 
also present for higher Reynolds numbers. Is the double threaded wake structure also 
present in the case of the highest Reynolds numbers, where the wake structure has 
a turbulent structure? If so, do instabilities in the wake cause kinking of the vortex 
threads which leads to this turbulent wake structure? Further research will address 
these questions in order to get a better understanding of the boundary layer separation 
from spheres at high Reynolds numbers as shown on this year's cover. 

4. Conclusions 

Flow visualizations of the wakes behind solid spheres moving under the action of 
gravity reveal remarkable differences with the wakes behind spheres held fixed: the 
crossing of threads of opposite-signed vorticity, the formation of kinks on these threads 
that develop into hairpin vortices. The ratio between the densities of the sphere and 
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(a) (b) (c) 

Figure 3. Stereoscopic views of rising zigzagging spheres and their wakes. The 
left part of each frame shows the xz-plane and the right part the yz-plane. The 
views illustrate the crossing at the center- line of the zigzag path of the two counter- 
rotating threads of the wake ('1'), the occurrence of kinks ('2') at the extremes 
of the path, and the formation of hairpin-like vortices ('3') as two neighbouring 
kinks connect. Values of the parameters d, ps/{p), G and Rej^ arc, respectively: 
(a) 5 mm, 0.95, 297, 450; (6) 5 mm, 0.95, 306, 475. As shown in (c) the crossing 
of the vortex threads results in a lift force L that is always directed towards the 
center-line of the zigzag path. D shows the direction of the drag force and B the 
one of the buoyancy. 




Figure 4. Sequence of stereoscopic views of a rising sphere and its wake. The 
left part of each frame shows the xz-plane and the right part the yz-plane. The 
views illustrate the process of formation of a hairpin-like vortex (d = 10 mm, 
p^/{p) = 0.99, G = 350 and Rbt = 576). 
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Figure 5. Phase diagram: density ratio Ps/{p) versus Galileo number. Tile grey 
box is the regime analysed by Jenny et al. [8] . They find in the most left region an 
axisymmetric wake. The symbols, directly taken from [8], denote: + steady and 
oblique, * oblique and oscillating regime with low frequency (0.045 < / < 0.068), 
X oblique and oscillating regime with high frequencies (/ = 0.180), O zigzagging 
periodic regime (0.023 < / < 0.035) and □ chaotic regime. The numbers denote 
the number of our experiment in table 1. Experiments 6 and 19 fall outside the 
diagram. 



that of the surrounding fluid appears to be important. Our experiments clearly show 
the difference in path and wake structure between rising and falling spheres with the 
same Galileo number. Furthermore, the double threaded wake structure seems to 
be a basic feature, even for large Reynolds numbers. This should be investigated 
thoroughly in future research. 
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